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Abstract 

We argue that it is possible to make a consistent picture of FNAL data including the 
production and decay of gluinos and squarks. The additional cross section is several 
pb, about the size of that for Standard Model (SM) top quark pair production. If the 
stop squark mass is small enough, about half of the top quarks decay to stop squarks, 
and the loss of SM top quark pair production rate is compensated by the supersymmet- 
ric processes. This behavior is consistent with the reported top quark decay rates in 
various modes and other aspects of the data, and suggests several other possible decay 
signatures. This picture can be tested easily with more data, perhaps even with the 
data in hand, and demonstrates the potential power of a hadron collider to determine 
supersymmetric parameters. It also has implications for the top mass measurement and 
the interpretation of the LEP Rb excess. 
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1 Introduction 



While there is still no compelling experimental evidence that nature is supersymmetric 
on the weak scale, there have been recent reports of data that encourage this view. The 
most explicit is an event in CDF m that does not have a probable SM interpretation, 
and can naturally be explained as selectron pair production^, |3|. In the interpretation 
when the lightest neutralino is the LSP, the analysis of this event leads to a fairly well 
determined range of masses and couplings for sleptons £, charginos Ci, and neutralinos iVj. 
Ci and Ni are the chargino and neutralino mass eigenstates, with C%, N% (C2, N4) being the 
lightest (heaviest) chargino and neutralino. If in addition there is a light stop squark t (i.e., 
m>i < mw)i it is remarkable that the chargino mass and couplings from || can explain 0] 
the LEP reported excess for Z — > bb decays (Rb), at least if that excess is not too large. 
Another encouraging result § is that the LSP resulting from these studies has a mass and 
coupling such that it is a good candidate for the cold dark matter of the universe, giving 
0.1 < n LSP /i 2 < 1. In the following, we use the Ni and Ci masses and couplings reported in 
Ref. §. 

If the stop squark t and at least one neutralino Ni are light, the top quark decay 
t — > tNi must occur along with t — » bW with branching ratio about one-half. Currently, 
the FNAL top quark pair production counting rate is interpreted as a measurement of 
a(pp — * tiX) x BR 2 (i — > bW). The reported values from CDF for this are about 7 pb [14], 
which is already larger than the predicted total cross section [J7j (about 5.5 pb) for the 
top quark mass (mt) extracted from kinematic reconstruction (about 175 GeV). The DO 
measurement for the production rate is smaller, but suffers from larger backgrounds. Unless 
there is additional production of top quarks or other particles with similar decay modes, 
there is no room for extra decay modes, such as t — * tNi, in this data sample (this is made 
quantitative below). 

In this paper we observe that super symmetry, with certain reasonable and well-motivated 
choices of sparticle masses, will lead to extra top quark production. Additionally, there are 
other final states from sparticles with a high purity of 6-quarks, leptons, and jets; they can 
mimic the top quark signal. If squark (q) and gluino (g) masses (excepting the stop squark) 
are between 200 and 300 GeV, then they have pb-level production rates at the Tevatron. 
We assume squarks are heavier than gluinos, motivated by the results of Ref. H as well 
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as the observation that gluinos decay predominantly to top quarks and stop squarks if the 
decays to other squarks are not kinematically allowed. A lower bound is set for the gluino 
mass by requiring m g > nit + n^i so that the decay g — > tt* can occur for physical masses. 
Since gluinos are Majorana particles, they will decay equally to tt* and tt ||. Pair produc- 
tion of gg will give a final state containing tt, tt, and tt in the ratio 1:1:2. If squarks 
are heavier than gluinos, they decay as q — > gq, Ciq, Niq; gq has the largest coupling but 
the smallest phase space. We assume that -R-parity is approximately conserved, as implied 
by the interpretation of the CDF event as sparticle production, so that all sparticles decay 
to Ni within the detector, and Ni subsequently escapes. The couplings to Cj and Ni are 
largely determined by the analysis of Ref. M. The size of the BR(g — > gq) is determined 
by the available phase space. We freely use the term Minimal Supersymmetric Standard 
Model (MSSM) to describe our models, and it should be understood that we do not assume 
gaugino mass unification, but we do assume that all squarks except the stop squark are 
mass degenerate. 

The relative masses of t and Cj are important for the signatures, and are not determined 
yet. Motivated by the LEP R\, data, we assume m t -<mc> 4 . Then the dominant decays 
are d — > bt*, and t — > cN ||. If half of the top quarks decay to stop squarks, then 
t{— ► bW)t{— » tNi) followed by t — > cNj, so half of all ti events give a Wbc signature; the 
finite detector acceptance, clustering of jets, etc., might lead to a similar final state in the 
Standard Model, but at a much lower rate. Similar observable differences occur in every 
distinct tt channel; results for some examples are shown in Sec. 4. Furthermore, a different 
value of nit is likely to be extracted from different channels if it is determined by comparison 
with a SM Monte Carlo, and the true value of nit might not be the apparent one. Note 
that a large number of charm jets arise from stop decays; if they could be tagged, e.g. by 
the lepton in charm semi-leptonic decays, it would help test our arguments. 

In the following section, we review several results on the top quark which are relevant to 
our discussion. In Section 3, we present a supersymmetric model which is consistent with 
the results of Section 2, explains several other pieces of data, and predicts several signatures 
in the present and future FNAL data samples. The detailed study of these models are shown 
in Section 4. Finally, Section 5 contains our conclusions. 
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2 Upper Bound on Non— SM Decay Modes of The Top Quark 

Based on a number of measurements and predictions, it is possible to bound non-SM 
decay modes of the top quark which are "invisible" to the standard searches. The final 
state resulting from tt — > bWcNiNi or ti — * cNiNicN\Ni would not have enough leptons 
or jets to be included in the leptonic, dileptonic, or hadronic event samples, and, hence, is 



invisible. Here, we review briefly a previous analysis bounding these "invisible" decays [10]. 

The FNAL experiments report essentially two independent measures of the top quark 
mass: the kinematic measure, whereby the four-vectors of all the decay products are re- 
constructed into the four-vector of the parent particle, and the counting measure, whereby 
one compares the observed number of events, corrected by efficiencies, to the production 
cross section as a function of mass. Additionally, based on single- and double-6-tagged 
events, CDF has reported a measurement of BR(i — > qW) for q ^ b. fll| . We interpret this 
as a limit on bw = BR(£ — > bW). The SM prediction for the production cross section at 
is bounded theoretically and the current best prediction of at at m t = 175 GeV is 5.52j/4 5 
pb from the first reference of 0. Given the measured quantities and the theoretical pre- 
dictions and their uncertainties, one can perform a x 2 minimization in the variables mt, at, 
and bw- Finding the minimum value Xmin yields mt = 168. GeV, a t j = 7.09l;gf pb and 
b w = 1.00l;5s- At the 95% confidence level, b w > -74, so an upper limit on BR(i — > X), 
where X ^ bW , exists of about 25%. 

This analysis disfavors a large component of non-SM top quark decays, such as a t — > 
cN%. However, this result has limitations. Namely, it does not include the possibility that the 
same physics which allows new top quark decays can also lead to more top quark production. 
(See [|l^] for an analysis which includes tt* production and the decays t — > bC\ and concludes 
that a light t is not excluded by the present FNAL data even for BR(i — > iiVj)=l/2). 
Furthermore, it assumes that all of the new decay modes result in final states which elude 
the standard searches. Finally, it does not include the LEP indirect fits to rat which favor 
a lower value (for example, the world average top mass including LEP, DIS, SLC, and 
FNAL data is 161±8 GeV assuming an 80 GeV Higgs boson, as expected from Ref. §). 
Since this would weight the x 2 f° r a smaller mt and a larger at, bw would be smaller, and 
more "invisible" decays would be allowed. Based on these observations, it is premature 
to conclude that a light stop is inconsistent with the observed top quark events. In the 
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following section, we present explicit supersymmetric models motivated by Ref. M which 
successfully address several of the loopholes in the previous analysis. Squarks and gluinos 
are produced at a significant rate and decay largely into top quarks. The stop squark is 
light enough that the charginos C\ and C<i decay to bt* , which in combination with other 
decays, give top-like final states. Finally, a lower top mass predicts a larger SM production 
rate, allowing more "invisible" decays. 

3 Top Quarks from Squark and Gluino Decay 

If selectrons, charginos, and neutralinos have masses of order mz, then squarks and 
gluinos might be light enough to be produced in significant numbers at FNAL. The analysis 
of Ref. H is done with a general low energy softly broken supersymmetry theory, without 
assumptions about gaugino or squark mass unification. As a result, the gluino and squark 
masses are not determined. However, there are phenomenological reasons to settle on the 
range 200-300 GeV for these masses. Given the analysis of Rb, we expect a light stop squark 
with mass less than about M]y. LEP and FNAL limits do not allow this to be too small [13], 



and we are forced to values somewhere between about Mz/2 and M\y. As explained 
in the previous section, light stop squarks alone dilute the signal t — > bW through decays 
t — > tNi at a level incompatible with the data, so a new top quark or top-like production 
mechanism is needed. The simplest method is to use the decay channel g — > tt* , which 
requires rrig > mt + m^. For this to be the dominant decay channel, the other squarks 
must be heavier than the gluino. Finally, the Tevatron has limited parton luminosities to 
produce heavy particles. Since about half of the top quark decays "disappear", we need 
a production mechanism which is about the same size as the SM rate. Therefore, we are 
naturally lead to squark and gluinos masses between 200 and 300 GeV. We will see that a 
number of observables depend on the particular masses, so eventually they can be directly 
measured. 

A similar mass hierarchy follows from theoretical considerations. Ref. ^ found gaugino 
masses obeying the mass relations M\ ~ M2, rather than the unification relation M\ ~ 
This could be explained by anomalous behavior of the U(l) mass, so that the non-Abelian 
masses may still approximately satisfy the unification relation M2 — M3. Then the gluino 
mass should be about three times the C\ mass, or in the range 195-270 GeV. Similarly, 
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we assume here that squarks are (except for stop mass eigenstates that are expected to be 
separated) approximately degenerate, and about 2.5 times the selectron mass, as suggested 
by models. For numerical work, we take a common squark mass nig for left- and right- 
handed squarks of five families which is slightly above nig. We study models with 160 < 
m t < 175 GeV, 210 < m g < 235 GeV, 220 < m 5 < 250 GeV, and 45 < m~ t < 60 GeV. 
All results are based on the analyses and models of ref. [§, ||] and are thus consistent with 
existing evidence for supersymmetry and with other particle physics constraints. 

4 Numerical Results 

In this section, we present separate results on the counting measurement of the top 
production cross section and kinematic measurement of the top quark mass in the super- 
symmetric models described in the previous section. All event simulation is performed using 
the Monte Carlo PYTHIA 5 . 7 with supersymmetric extensions |l5| . The cross sections are 
computed at the Born level with additional QCD radiation added in the leading-log ap- 
proximation. The structure functions used are CTEQ2L, and the hard-scattering scale 
used in the evaluation of the structure functions and the running couplings is the partonic 
center of mass energy (transverse mass) if the partonic final state has no (some) non-zero 
QCD quantum numbers. Particle energies are smeared using Gaussian resolutions based 
on the CDF detector, and jets are defined using the PYTHIA LUCELL subroutine. Jets are 
6-tagged with a constant efficiency when they contain a high-py 6-parton; the exact effi- 
ciency is stated only when a result depends upon the choice. All leptons and photons must 
be isolated from excess transverse energy. 

To eliminate dependence on the particulars of 6-tagging and isolation efficiencies and 
detector cracks, we will present some results as ratios with the SM signal expected using 
the same Monte Carlo routines. Where there is no SM signal expected, or just a small 
one, then we show an expected number of events in 100 pb _1 . First, we examine the 
counting measurement of the top quark production cross section. Since there are so many 
production processes and decay chains at work when all sparticles have weak scale masses, 
it is clear that a multi-channel analysis should be performed. To illustrate this point, we 
present a Unitarity Table (Table 1), based on a very loose set of experimental cuts, so that 
almost 100% of the simulated events fall into some category. On one hand, this is useful 
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to check that a particular model does not contradict data by predicting an anomalously- 
high rate in some channel. On the other hand, it can point the way to unexpected features 
in the data, or explain why various signals are not present. This particular table is for a 
representative model, but displays the essential features. The columns represent various 
production mechanisms: (1) SM tt, (2) ti with supersymmetric decays, (3) squark pair, (4) 
squark-gluino, (5) gluino pair, and (6) d or Ni in association with a squark or gluino. The 
rows represent various final states defined by the presence of charged leptons photons 
7, and a large missing transverse energy $t) and the number of jets rij (if present, 
and 7 are explicitly noted) : (a) rij < 3, (b) rij > 3, with small (c) rij > 3, with 
large Ifrr, (d) the subset of (b) and (c) which has rij > 6, (e) ^,-rij < 3, a lepton with less 
than 3 jets, (f) f^.rij > 3 and large $ T , (g) ^.nj > 3 and smaller, (h) ^±,7^- > 
with large IpT, (i) (^7, rij > 0, (j) 7, rij > 0, and (k) 77, rij > 0. The numbers in each 
column represent the fraction of generated events in a particular final state. If the most 
important final states are accounted for, then the sum for each column should be close to 
unity, excluding the 6 jet final state, which is a subset of two others. Indeed, the sums vary 
between .98 and 1.00. This table demonstrates that we understand where the individual 
supersymmetric contributions go. Only after multiplying by the production cross section 
(and choosing more realistic cuts) for each process can we determine the observable rate. 
Row (c) corresponds to the standard SUSY search mode of multi-jets plus $t- For the 
particular choices of squark and gluino masses considered here, we still elude the present 
experimental bound. 

In the following, we present results based on more realistic cuts. There are three standard 
top quark search modes defined by the CDF cuts: (i) leptonic, (ii) dileptonic, and (hi) 
hadronic. We find substantial signals in 2 additional channels: (iv) "W"bc, and (v) 76c. 
The channel (iv) cuts are the same as for (i), except only 2 jets are allowed. The excess 
number of " W" be events (where " W" may have a different transverse mass since Ni 's carry 
away energy) would appear as an excess of W plus 2 jet events with one 6-tag; the second 
jet is charm which can be tagged with a lower efficiency. The channel (v) cuts require a 
high-pt 7 (> 20 GeV) in the central rapidity region (|r/ 7 | < 1), with 2 or more additional 
jets. One of the two leading jets must have a 6-tag. As discussed earlier, channel (iv) should 
be limited in the Standard Model. Channel (v) events arise from the decay N2 — ► iVi7 in 
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Table 1: Unitarity Table, illustrating the fraction of events which fall into several categories 
for SM ti and MSSM processes. 



association with top quark or top-like decays, and should have a tiny contribution from ti 
production alone in the Standard Model. 

Table 2 summarizes the results of the counting experiments for the various models 
considered. The numbers have ranges, and we comment on the correlations later. In the 
upper portion of the table, we present fractions in each row of our Monte Carlo estimate of 
the process rate after cuts divided by the same estimate for SM top quark production. These 
numbers suggest that different values will be obtained for the top production cross section 
in different modes; in a given mode, the value will depend on the analysis and cuts. This is 
consistent with the reported CDF and DO cross sections [14]. The row labelled £^£^ shows 
the predicted non-SM signal of like sign leptons; this is possible because of the Majorana 
nature of the gluino. We note that about 1/7-1/5 of all dilepton events should have leptons 
with the same charge. The middle section of the table shows the expected number of events 
in 100 pb _1 for the two aforementioned channels "W"bc and jbj. These numbers do not 
include a 6-tagging efficiency Also included in the 'ybj sample is the expected number of 
events from Cj(— > bt)N2(^ N\j) production (+35). The final section shows the variation in 
total production cross sections for the various channels. Note that the MSSM ti production 
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Table 2: Expected results of the top quark counting experiments for the MSSM. The ap- 
parent top production cross sections are shown in the final column. The number of events 
in the present data sample for two channels are displayed in the middle section. Typical 
MSSM production cross sections appear in the final section. 



cross section is identical to the SM one; the only difference occurs in the allowed top quark 
decays. Of course, the various apparent cross sections are correlated. For m t =160 GeV, 
mq=220 GeV, m g =21Q GeV, m x o=38 GeV, m^=45 GeV, the cross sections measured in 
the three modes are 6.5, 7.8, and 6.6 pb. For m t =165 GeV, m,j=240 GeV, m,g=220 GeV, 



in 



,o=38 GeV, 771^=50 GeV, the numbers are 5.7, 5.3, and 6.3 pb. 



Some of the larger apparent rate for top quark production for SUSY processes comes 
from the increased cross section for smaller mj. We use the resummed prediction, which is 
similar to the NLO number. The production of squarks and gluinos, on the other hand, are 



calculated only at LO, and could receive a substantial NLO correction. Based on Ref. [16], 
we estimate that our squark and gluino production cross sections (which are evaluated using 
the transverse mass as the factorization scale) could be increased (but not decreased) by as 
much as 40%. We have not included this ivT-factor in our study. In addition, smaller mt 
allows for smaller gluino and squark masses, which further increases the MSSM rate. We 
have made no attempt to optimize the numbers in Table 2. It is remarkable how naturally 
the apparent cross section values span the experimentally allowed values. 
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Figure 1: The distribution of "ybj events expected in 100 pb . There are contributions 
from the CiN 2 , tt, and q,g production processes. 

Since they are a novel feature of our models, we also show a typical scatter plot of the 
events expected in 100 pb -1 with signature 7&$r+jets in Fig. 1, resulting from models with 
m t =160 GeV. They are of particular interest because there is no parton-level SM source 
of such events, and our MSSM scenario predicts a significant number. There are three 
sources: qq Ci(-» bt*)N 2 (^ -yN^qq^gg) -> t(-> bW(-+ t*N 2 (^ 7JV1)), with 

t — > cJVi, and cascade decays from qq, qg, gg, gN,i,gCi, q~Ni,qCi, populating different regions 
of the plot. The C{N 2 and tt signals depend mostly on the super symmetric interpretation 
of the CDF event and the postulate of a light stop squark to explain R^; their signal may 
be present regardless of the other squark and gluino masses. Note that the first of these 
produces only two prompt jets, while the other two produce several jets. Finding these 
events could confirm supersymmetry in general and our arguments in particular. 

In addition to the top quark measurements based on counting events, there are kinematic 
measurements, such as the reconstructed top mass, the tt invariant mass, the transverse 
momentum of the pair, etc. The explicit reconstruction is a difficult task. Not all of the 



10 




Figure 2: The invariant mass distribution of the tt pair for SM and MSSM top quark 
production compared with the CDF data. 



11 



apparent top quark signal in our models comes from real top quarks, but a substantial 
fraction does. Rather than attempt an explicit reconstruction for such detector dependent 
quantities, we have tried to determine where we should expect agreement with and devia- 
tion from the SM distributions. We identify the partonic top and anti-top quarks in our 
generated events and use these to calculate the invariant mass and transverse momentum 
of the ti pair. In the process, we ignore those events which have only one or no real top 
quarks, but scale those events with two top quarks to the full rate including the discarded 
events. The invariant mass distribution from SM top quark production alone with m^=170 
GeV, for our MSSM with mi=160 GeV, and with the CDF data is displayed in Fig. 2. The 
CDF data is from [|17]]. While we have only performed a crude simulation, we believe this 
figure demonstrates that our MSSM is consistent with the data. This consistency is under- 
standable. Since top quarks are coming from gluino decay just above threshold, they are 
produced almost at rest in the lab frame. As a result, the distribution must peak slightly 
above 2mt- 

The transverse momentum of the ti pair is displayed in Fig. 3 for the same conditions 
as for Fig. 2. Here, we have not smeared the distributions, but we observe that the MSSM 
distribution is broader than that expected in the SM. This too is expected; the two-body 
decay of the gluinos tends to randomize the top quark direction, removing the approximate 
balance in pt expected from SM ti production. The expected SM ti distribution displayed 
in Ref. jL7j] is narrower than the data, so the MSSM could explain this discrepancy. 

The events used for the kinematic reconstruction of mj with a 6-tag come from the 
^z/)+jets mode. The lepton and the IpT hi these events should have a transverse 
mass consistent with that from the decay of a W boson. When top quarks are produced in 
MSSM events, there can be additional $r from cascade decays down to N\. The expected 
distributions for the transverse mass and the CDF data are shown in Fig. 4. The naive 
expectation that the MSSM distribution must be distinguishably different is not fulfilled. 



The CDF data is from |8|. 

Since many of the top events have associated jets, the apparent top mass deduced 
from such events will only be the actual top mass if very particular cuts and analyses are 
used. For example, we have noticed that the invariant mass distribution of the leptons in 
dilepton events is softer than for SM top events. This indicates that the mass kinematically 
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Figure 3: The transverse momentum distribution of the tt pair for the SM and MSSM 
compared to the data. The generated distributions have not been smeared. 
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Figure 4: The transverse mass of the lepton and $t for the SM and MSSM compared to 
the data. 
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reconstructed from dilepton events will be lower in the MSSM than for the other modes. 
Note that this is the only mode which does not require a 6-tag, so that the additional jets 
from squark decays can enhance the signal. 

5 Conclusions 

We have argued that existing data is consistent with the possibility that hundreds of 
squarks and gluinos have been produced at FNAL. Squarks decay mainly into gluinos, 
charginos, and neutralinos, gluinos into top quarks and stop squarks, and BR(i — ► tiVj) is 
about 1/2. We have checked that the predicted counting measures and kinematic measures 
are consistent with the available data, and, in some cases, give a better description. A 
number of associated predictions allow this view to be tested, possibly with existing data. 
If correct, it has implications for the top quark mass and cross section measurements, for 
interpreting the LEP i?5 data, and of course for the existence of supersymmetry in nature. 
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